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ABSTRACT 

The integration of embedded generation (EG) units into the distribution system brings technical and economic 
benefits to the power system network. But, placement of EG's at non optimized locations can effect voltage profile with 
increased system losses. This paper intends the optimal placement and contingency analysis of several Embedded 
Generation (EG) units in distribution network using Genetic Algorithm (GA). The main intention is to diminish the total 
real power losses with proper voltage profile. It go behinds the conversation of evaluating the effect of the location and the 
size of EG to the system before and after contingency is created in the system by fault. The projected allocation method 
and contingency analysis is verified using IEEE-33 and 69 bus radial distribution systems. The algorithm has been 
programmed in MATLAB and results are compared. 
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INTRODUCTION 

In recent years, due to the increasing interest on renewable sources such as wind, solar, geothermal, 
biomass and ocean energy, all over the world, the number of studies on embedded generation (EG) has rapidly increased 
and makes EG units cost effective. By this, injection of embedded generation systems to distribution networks not only 
diminishes the power losses but also improves voltage profile of the network. 

The problem of global warming and exhaustion of fossil fuel have fostered the increased usage of less 
environmentally-polluting embedded generation (EG) sources closer to load demands in the distribution system [5]. 
EG units are rapidly increasing and becoming an impressive alternative to grid for reliability and ability to exploit 
renewable energy sources. 

Literature survey shows that a lot of work has been done on optimal location and size of EG units in distribution 
networks. Some of such studies are: Tabu search [1], Artificial Intelligence (AI), Evolutionary computation, 
fuzzy-GA algorithm, Ant Colony Search Algorithm (ACS A) [2], Particle Swarm Optimization (PSO) are some of popular 
optimization techniques that are normally used to solve this type of problems. The associated works is carried on GA for 
optimal placement of EG using MATLAB as a tool for solving that problem which has been proposed in. In additionally 
we worked on multiple EGs with and without contingency analysis. 

This paper provides solution for problem of optimal allocation and contingency analysis by using Topology based 
load flow solution technique [3] with genetic algorithm by considering EG as PQ model. When compared to other search 
techniques in this proposed technique convergence is fully guaranteed and computation time is very Minimum 
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TOPOLOGY BASED LOAD FLOW SOLUTION 

In Distribution networks R/X ratio is very high and due to this, these networks are unconditional and conventional 
therefore Newton-Raphson (NR) and fast decoupled load flow (FDLF) methods are inaccurate in providing solution for 
such type of networks. Here in this paper, Topology based approach is used [4]. This proposed technique contains only 
evaluation of a simple algebraic expression of voltage magnitude without caring of trigonometric terms so that it eliminates 
the complex process of identifying nodes connected beyond a particular node. The two matrices, bus injection to node 
power matrix [BINP] and line loss to node power matrix [LLNP] are created based on connectivity of network. Under 
assumption that shunt capacitance is negligible at the distribution voltage level chance of convergence is fully possible. 

It is assumed that the three-phase radial distribution networks are balanced. This assumption is quite valid for 
1 lkV rural distribution feeders in India and elsewhere. 

Formation of BINP, LLNP Matrices 

Let us Consider 6- Bus Radial Distribution Network (figure 1) 
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Figure 1: Sample 6-Bus Radial Distribution System 
Bus-Injection to Node Power Matrix (BINP) Matrix 

For fig 1.1.1, N2, N3, N4, N5 and N6 are the equivalent powers at each Node. 
N2= S3+ S4 + S5 + S6 
N3 = S3 + S5 
N4 = S4+ S6 

Where S3, S4, S5 and S6 are complex load powers at buses 3, 4, 5 and 6 respectively 

[N] = [BINP] [S] (1) 
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Line Loss to Node Power matrix (LLNP) matrix: The relation between the line losses and node power can be 
obtained by following equations: 
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N2 = SL1 + SL2 +SL3+ SL4 +SL5 
N3' = SL4 
N4'= SL5 

Complex power losses associated with lines 1, 2,3,4,5 are given by SL1, SL2, SL3, SL4 and SL5 respectively. 
Sum of the line losses appearing at Node 2, 3, 4, 5 and 6 are given by N2 , N3\ N4\ N5\ N6\ 



[NT] = [LLNP] [SL] 
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Effective load at each node = N + N (3) 
The algorithm to solve the load low problem is given below. 

• Read the system data, Vl= 1.0 p.u. Line losses are assumed to be zero in the first iteration. 

• Build BINP matrix and LLNP matrix. 

• Obtain P effect i V e + j *Q effective, at each node. 

• Initialize iteration count =1 

• Obtain receiving end voltages using V2=V (B[j]-A[j]) and power loss on all lines using PlossUI 
=RU]*(P 2 2 +Q2 2 )/V 2 2 , QlossD] =XLj]*(P 2 2 +Q 2 2 )/V 2 2 and S L0S s =Ploss + JQloss 

• Multiply the power loss column matrix, SL with LLNP matrix to get N' matrix. N' represents the part of line 
losses in the effective load at various nodes. [N'] = [LLNP] [SL] 

• Calculate the total effective load at various nodes by adding the N and N' matrix. Effective load at each node = N 
+ N' 

• Increment the iteration count Repeat the steps from step (5) using new effective loads at every node. 

• If the difference in the voltages between present iteration and previous iteration is greater than 0.001 pu, then 
increment the iteration count and repeat from step (5), otherwise, print the result. 

FUNCTION OF GA IN OPTIMAL PLACEMENT AND SIZING OF EG 

The genetic algorithm (GA) is an optimization and searching technique which operates on the principle of 
genetics and natural selection [6], [7]. Here, continuous or floating GA is used as it is having a benefit of exact 
representation of the continuous parameter. 

Continuous Genetic Algorithm (GA) is used to determine the optimal placement and size of Embedded 
Generation units in the distribution system in order to diminish the total losses in the system. By reducing these losses, the 
voltage profile at each bus is expected to be enhanced [8]. This scheme needs load flow to run for several times. Later, 
getting the optimal placement and size simultaneously for Embedded Generation (EG), the algorithm is finished. 
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The objective function is the results of total loss of the system, Ploss to be minimized, H as follows: 

Where, nline is number of lines in the system, Ploss is the real power losses at each bus, H is the total real power 
losses that are to be minimized. 

T he total length of chromosome is equal to (number of EG's x2). the first variable represents the location 

and the second one represents the size of EG in chromosome. 
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Figure 2: Chromosome 

After evaluating each chromosome, the objective function is transformed and normalized to a fitness scheme to be 
maximized as follows: 

f=l/ (1+H), f is the fitness function that is to be maximized, v 

CONTINGENCY ANALYSIS 

Contingency analysis is the study of the outage of elements such as transmission lines, transformers and 
generators, and investigation of the resulting effects on line power flows and bus voltages of the remaining system. It is 
important for both the system operator and the system planner to be able to evaluate how the line flows and bus voltages 
will be altered in the new steady state. Overloads due to excessive line currents must be avoided and voltages that are too 
high or too low are not acceptable because they render the system more vulnerable to follow on cascading outages. The 
large numbers of possible outages are analyzed by means of contingency analysis. Great precision is not required in 
contingency analysis since the primary interest knows whether or not an insecure or vulnerable condition exists in the 
steady state following any of the outages [9] . 

This chapter focuses on the IEEE-33 and IEEE 69-bus test systems on which the analysis is done by introducing 
the contingency situation which is referring to the fault (i.e. line 3- 28 in IEEE-69 and line 2-19 in IEEE-33) to see the 
impact of EG that has or have been installed in the system. The study will emphasize on the changes of voltage profile and 
the system losses before and after reconfigurations caused by the fault with EGs. However, the type of faults is not to be 
considered since it is assumed that the system is reconfigured after the isolation of the fault. This study will help in finding 
the trends of optimal size and location of EG's in these test systems. 

OPTIMAL ALLOCATION STUDY AND DISCUSSIONS 

The proposed method has been tested on IEEE-33 and 69 bus radial distribution systems. The algorithm has been 
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programmed in MATLAB and the test system can be obtained. 
The GA properties are as follows: 

• Mutation probability, pm = 0.5 

• Selection=0.5 

• Population = 50 

• Number of EG units = up to5 

The reason for using 5 units of EG's that there is no considerable reduce in real losses when adding EG units more 

than 5. 

Variation of Total Power Losses with Respect to Iterations for IEEE-33 and 69 Bus Systems without Contingency: 

From the figure 3 it is clear that, for IEEE-33 bus system the Real power losses without placing EG is 0.308PU, 
while connecting 5 EG units, active power losses are reduced to 0.0195 PU. The total real power losses are decreased by 
93.66% and for IEEE-69 bus system minimum active power losses of the given distribution network are 0.0205PU while 
connecting 5 EG units at optimal locations the total real power losses are reduced 0.0063PU. 
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Figure 3: Total Real Power Losses Vs Iterations for IEEE- 33 and 69 
Bus Systems with EG's and without Contingencies 

Table 1: Effect of Multiple Placement of EG Units without Contingency on IEEE- 33 and 69 Bus Systems 



No of EG 
Units 


IEEE 33 Bus System 


IEEE 69 Bus System 


Embedded Generation 


Total Real 

Power 
Losses(PU) 


Embedded Generation 


Total Real 

Power 
Losses(PU) 


Size (PU) 


Locations 


Size (PU) 


Locations 


No EG 






0.308 






0.206 


1 EG 


1.5338 


13 


0.1647 


0.6464 


63 


0.0215 


2 EG 


1.0936, 1.0846 


19,30 


0.0707 


0.5672,1.1962 


69,63 


0.0162 


3 EG 


1.1230,1.7204, 
1.1465 


9,19,20 


0.0453 


0.4121,1.1132,0.0772 


61,6,47 


0.0107 


4 EG 


1.0921,0.9882,1 
.1023,1.0831 


10,13,16,1 
8 


0.0259 


1.4898,1.3809,0.4098, 
1.2123 


59,64,47,5 
8 


0.0084 


5EG 


1.0696,0.9949,1 
.1127,0.9585,1. 
3856 


10,18,22,2 
6,31 


0.0195 


0.8143,0.5173,0.5882, 
0.1535,0.6718 


38,49,51,5 
3 ,63 


0.0063 
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Voltages at Different Buses for IEEE- 33 and 69 Bus Systems without Contingency 

From figure 4 It can be observed that when EG units are inserted in distribution network the voltage profile is 
improved from without EG to with 5EG units. 
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Figure 4: PU Voltages at Different Busses for IEEE- 33 and 69 Bus Systems with and Without EG's 

Variation of Total Power Losses with Respect to Iterations for IEEE- 33 and 69-Bus Systems with Contingency 

From the figure 5 it is clear that ,for IEEE-33 bus system, the Real power losses without placing EG is 0.2606PU, 
while connecting 5 EG units, active power losses are reduced to 0.0196 PU. The total real power losses are decreased by 
98.2% and for IEEE-69 bus system, minimum active power losses of the given distribution network are 0.196PU while 
connecting EG units at optimal locations the active power losses are reduced to 0.0156PU. 




Figure 5: Total Real Power Losses Vs Iterations for Ieee-33 and 69 Bus 
Systems with EG's and with Contingencies 



Table 2: Effect of Multiple Placements of EG Units with Contingency on IEEE- 33 and 69 Bus Systems 



No of 
EG 
Units 


IEEE 33 Bus System 


IEEE 69 Bus System 


Embedded Generation 


Total Real 

Power 
Losses(PU) 


Embedded Generation 


Total Real 

Power 
Losses(PU) 


Size (PU) 


Locations 


Size (PU) 


Locations 


No EG 






0.2606 






0.196 


1 EG 


1.2953 


11 


0.1362 


1.2537 


61 


0.0426 


2 EG 


0.9687,1.2604 


16,27 


0.066 


1.5237,0.5326 


61,39 


0.0298 
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Table 2: Contd., 


3 EG 


0.7429,1.6654,1.1225 


18,25,27 


0.0424 


0.2352,1.4387, 
0.1133 


4,27,50 


0.0203 


4 EG 


0.5665,1.3878, 
0.6069,1.6769 


8,14,16,25 


0.0231 


0.0320,1.2650, 
1.1938,0.5099 


3,15,50,51 


0.0194 


5EG 


1.7000,1.5801, 
0.465,0.8817,1.0706 


6,14,16,21,28 


0.0196 


1.4591,0.5851, 
0.5798,0.559,0.5672 


3,4,56,61,24 


0.0156 



Voltages at different buses for IEEE 33 and 69 bus systems with contingency: 

From figure 6. It can be clearly observed that when EG units are inserted in distribution network the voltage 
profile is improved from without EG to with 5EG units. 




BUS NUMBER BUS NUMBER 

Figure 6: PU Voltages at Different Busses for IEEE- 33 and 69 Bus Systems with Contingencies 

From the table 3, it is clear that the real power losses are reducing when EG is placed at optimal location and with 
optimal size. The percentage improvement from without EG to 5 EG is around at 90 in every case i.e. with contingency 
and without contingency. 

Table 3: Comparison Table 



Type of 
Standard 
Bus 




Total real power losses (PU) 


Without 
EG 


5EG 


% 

Improvement 


IEEE-33 

Radial 

Bus 

System 


Without 
contingency 


0.308 


0.0195 


93.66 


With 

contingency 


0.2606 


0.0196 


92.47 


IEEE-69 
Radial 
Bus 
System 


Without 
contingency 


0.2063 


0.0063 


96.94 


With 

contingency 


0.1960 


0.0156 


92.04 



CONCLUSIONS 

This paper provides a continuous GA for determining optimal placement and sizing of EG units. Here, Multiple 
EG installation cases were studied using IEEE 33 bus distribution system. The results demonstrated and emphasized that 
multiple EG installations decreased total real power losses more than single EG installations. However, multiple DG 
installations may result in unnecessary additional costs. In this paper "Topology based approach for Load flow Solution of 
Radial Distribution Networks" is used which is simple, efficient and accurate method. 
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